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1. Introduction 

Since the first observation of Bose-Einstein condensation (BEC) in trapped dilute 
atomic gases [1112], developments of experimental and theoretical studies in this field 
have advanced rapidly |S1 Ej- One of the greatest advantages of these systems is 
that experimentalists have useful tools to create various external potentials for atoms. 
Magnetic traps, which are standard equipments to confine atomic gases, create harmonic 
potentials for atoms. In early stages of the study, many characters of BECs in harmonic 
confinement were vigorously investigated. The laser beams can also make periodic 
potentials for atoms, called optical lattices, by superposing the counter-propagating 
laser beams and preparing standing waves j3]. Moreover, the blue-detuned laser beams 
act as potential barrier for atoms, and one can create double-well potentials by focusing 
the laser beam into the center of a confining potential One can regard Bose-Einstein 
condensates in a double-well potential as a bosonic superfiuid-insulator-superfluid (S-I- 
S) junction like a fermionic superconductor-insulator-superconductor junction. 

The controllability of external potentials provides a suitable stage for study 
of quantum tunneling phenomena. Actually, many theorists have investigated the 
Josephson-like effect, corresponding to a kind of quantum tunneling |31 IH Ej, since 
the experimental realization of BECs in a double- well trap jH] . Recently, some types of 
the Josephson effect, called Josephson plasma oscillation and self-trapping in Ref. [Zj, 
were experimentally observed [S]. 

The study of elementary excitations is one of the main subjects to understand 
various properties of the Bose condensate, such as dynamics, thermodynamics and 
superfluidity. Several interesting properties of the elementary excitations have been 
probed experimentally in gaseous BECs, including the observation of collective modes |31 
and the measurement of the Bogoliubov excitation spectrum with use of Bragg 
spectroscopy |3l El ^j- As for the tunneling problem of the elementary excitations, 
it was predicted that low-energy excitations perfectly transmit through a potential 
barrier, and such a behavior was called the anomalous tunneling jTHl [Ej. A recent 
work shows that the anomalous tunneling property is crucial to the phononlike form of 
the excitation spectrum of BECs in a periodic potential, which is directly connected to 
the superfluidity [TK] . 

On the other hand, the macroscopic phase difference ip plays a key role in Bose 
condensed systems. Experiments on the matter wave interference jHj and the Josephson 
effect [S] have clearly demonstrated the existence of the phase difference between two 
spatially-separated BECs. In addition, the problem of the tunneling of quasiparticles 
in the fermionic superconductor-insulator-superconductor junction has been studied in 
detail, and it is well known that the phase difference between two superconductors 
essentially affects the scattering processes at the potential barrier . It is also expected 
that the tunneling of the elementary excitations in bosonic S-I-S junctions depends on 
the phase difference. 

In the present paper, we study the tunneling of the elementary excitations in a 
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bosonic S-I-S junction and analyze the <y9-dependence of the transmission coefficient T 
by solving the Bogoliubov equations analytically. It is found that T strongly depends 
on If, with regard to the anomalous tunneling. The peak width of the transmission 
coefficient decreases as (p increases, and the peak vanishes when the phase difference 
reaches — f which gives the maximum supercurrent of the condensate. We will 
show that the anomalous tunneling and its remarkable yj-dependence originate from the 
existence of the components localized near the potential barrier. 

The outline of the paper is as follows. In Sec. II, we introduce a formulation of the 
problem using the mean field theory. We calculate the condensate wave function and 
the relation between the supercurrent and the phase difference analytically. In Sec. Ill, 
we analytically solve the Bogoliubov equations and obtain the transmission coefficient. 
In Sec. IV, we discuss a mechanism of the anomalous tunneling. We also discuss the 
feasibility of observing the anomalous tunneling in real systems. We summarize our 
results in Sec. V. 



2. Mean field theory and model potential 

We consider a EEC at the absolute zero of temperature in a box-shaped trap which 
consists of a radial harmonic confinement and end caps in the axial direction. We assume 
that the frequency uj± of the radial harmonic potential is large enough compared to the 
excitation energy for the axial direction. Then, one can justify the one- dimensional 
treatment of the problem. Such a configuration was realized in a recent experiment |17j . 
It is assumed that the axial size L of the system is so large that the effect of the edge 
of the system can be neglected. Setting a potential barrier at the center of the BEC, 
one can create a bosonic S-I-S junction. In order to treat the problem analytically, we 
adopt a 5-function potential barrier as the potential barrier, 

V{x) = Vo6{x). (1) 

The schematic picture of the bosonic S-I-S junction is shown in Fig. lAll 

Our formulation of the problem is based on the mean field theory, which consists 
of time-independent Gross-Pitaevskii equation and the Bogoliubov equations. They are 
written in the dimensionless form as 



2dx 

and 



[-1^, -U + V{x) + \^o{m^o{x) = 0, (2) 



— _ 2 V / \ / 

Hq —"^q^x) \ I U{x) \ _ [ u{x) 



\l'o(a;)*^ —Hf) J \ v{x) J \ v{x) 



(3) 



Ho = - ^ + V (x) + 2|'l'or. (4) 

We have introduced the following notation: 

^ = 7> (5) 
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i^o,u,v) = J — {^o,u,v), (6) 
V ^0 



/i = ^, (8) 



V{x) = Vo6{x) = ^6{^), (9) 

where fi is the chemical potential, Uq is the density of the condensate fraction for x ^ ^, 
and e is the excitation energy. The healing length ^ is expressed as = --^^=. Here, 
\E'o(a;) is the wave function of the condensate and {u{x),v{x)y is the wave function of 
the excitation. The coupling constant g is affected by the harmonic oscillator length a± 
of the radial confinement as g = "^^^ [TH] . We shall omit the bars for all variables in 
Eqs. ©-(jni) hereafter. It is clear in Eq. (jH)) that the density of the condensates acts as 
a kind of potential for excitations. 

Since the purpose of our study is to investigate the relation between phase difference 
and the tunneling of elementary excitations, we need to obtain the condensate wave 
function with phase difference from Eq. It corresponds to a solution of Eq. (j21) 

which has stationary Josephson current. Let us first find such a solution assuming 

^o(x) = e'(^"+^±) (10) 

at far from the potential barrier |x| ^ 1. The constant C± expresses the phase of the 
condensate wave function at x ^ ±oo. If there is no interatomic interaction, or (7 = 0, 
the solution of Eq. Q satisfying the boundary condition (fTUI) does not exist. This 
is because finite fraction of the incident wave is inevitably reflected by the potential 
barrier. On the other hand, if there is repulsive interatomic interaction, the solution 
with the boundary condition (fTUI) exists as shown in Refs. fOl 120] ■ That means that 
BECs with repulsive interaction go through the potential barrier without reflection, and 
this behavior clearly exhibits the superfluidity of the BECs. In this case, the chemical 
potential is expressed as 

/i = i + Y- (11) 

One can realize such a situation in experiments by moving the potential barrier at the 
velocity of — ^ and choosing the coordinate system with their origin at the potential 
barrier. Substituting = A{x)e^^^^\ we can rewrite Eq. (j2)) as 

- i (0 - I'A-'^ + (V(^) -1-'^)A + A' = 0, (12) 

where q expresses the current of the condensate fraction in dimensionless form and 
Eq. ()13p corresponds to the equation of continuity. In dimensional form, the current 
is expressed as Ml!o_ We do not explicitly write V{x) in Eq. (fT^ hereafter, because 
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the effect of the 5-function potential barrier appears only in the boundary condition at 
X = 0. Multiplying Eq. (fT^ by jA, one can integrate the equation as 



1 
4 



dx 



_ ^A'^ - (- + ^)A^ + -A' 
4 4 4^ 4 



C, 



(14) 



where C is an integration constant. The integration constant can be determined by the 
boundary conditions 

(15) 
(16) 



A{±oo) 
dA 



dx 



itoo 



1, 

0. 



Then, Eq. fll4|) can be written as 



A 



dA\- 
dx J 



1 - A'YiA' - q' 



Integrating Eq. (fTTj) again, one obtains 



where 



A' = ^(xf + 



7(x) = \/l — g^tanh ( y'l — g^(|x| + xq 



[17) 



:i8) 



(19) 



Substituting Eq. (fTH|) into Eq. (fT^ . one can obtain the phase S{x) of the condensate 
wave function, 



S{x) - S{0) 



qx + sgn(x) 



tan^ 



^1 / li^y 



tan 



^J7(0)' 



(20) 



One can easily find from Eqs. (fT5|) and ((201) ^^lat the condensate wave function is 
expressed as 



^q(x) = ei(9^-«g'^(^)''o) (7(0;) - sgn(x) iq) 



where 



7(0) - ig 
7(0)2 + q2 



(21) 



(22) 



This solution is almost the same as that for a gray soliton [21], which is a kind of 
nonlinear excitations of BECs with repulsive interaction. The only difference from the 
solution for a gray soliton is the constant Xq in 7(0;) which depends on the potential 
strength Vq. The constant xq can be determined by the boundary conditions at x = 0, 

^o(+0) = *o(-0). 



d'^r 



dx 



+0 



(23) 

.o + 2Fo^o(0). (24) 
Substituting Eq. (PT|) into Eq. (j211), one obtains the equation to determine Xq, 

7(0)^ + Vo7(0)' - (1 - g')7(0) + Voq' = 0. (25) 



dx 



Phase Dependence of Phonon Tunneling in Bosonic S-I-S Junctions 



6 



Solving Eq. (j25p . one obtains function of q and Vq. The solution (j2H) has been 

obtained in Ref. |19j . 

We shall next define the phase difference (p and express the current q and the 
constant Xq (or 7(0)) by using ip. According to Ref. j22I, the phase difference is given 
by 



p^q I dx{ — -l). (26) 



Substituting Eq. into Eq. (j^ . one obtains 



V? = 2 



tan-f^Vtan-f2M' 



(27) 



\ (1 / \ (1 , 

Equations and (j?fj) yield q and 7(0) as functions of Vq and p. Assuming Vq ^ 1, 
we expand Eqs. and (jTfj) into power series of and obtain 

simp f cosp 2 + COSV5 — 3cos^(y9\ , , 

^(0) ~ ^ + ^Q^'^ (i _ ^-^Q^^ _ l+4cosy.-3cosV \ 

The leading term of Eq. (j^Hj) is the well-known relation between the Josephson current 
and the phase difference. It is obvious from Eq. that there is the critical current, 

2k ~ 

at the critical phase difference, 

TT 1 1 

Pc^ + — ^- (31) 

As the potential strength increases, the critical current decreases. The leading term in 
Eq. (jHUIl is consistent with the result in Ref. [12]. We can easily see from Eqs. and 
()30|) that the critical current q^ equals to zero when (7 = 0. 



3. Calculations and results 



In this section, we shall solve the Bogoliubov equations with the condensate wave 
function of Eq. (j2H) which has the phase difference and discuss the tunneling of the 
elementary excitations. 



3.1. Definition of transmission and reflection coefficients 

We need to clarify the general definition of transmission and reflection coefficients, before 
we discuss the tunneling of the elementary excitations. In the case of the tunneling of a 
single particle, it is well-known that one can easily define the transmission and reflection 
coefficients by means of the equation of continuity for the probability current. When 
one discusses the tunneling of a single particle through a potential barrier l^b(a;), one 
usually considers a situation in which the particle comes from the left (or right). In this 
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situation, one assumes that the solution 4>sp{x) of the Schrodinger equation takes the 
form of 



M^)=\ ^ . ,L ' (32) 



X 



where the wave number k± is related to the energy of the particle E as 



2m 

fc±^^^(E-ypb(±oo)). (33) 

The equation of continuity for the probability current gives the relation: 

|aspP + ^|&spP = l. (34) 

This equation means that the transmission coefficient Tgp and the reflection coefficient 
i?sp are defined as 

k+ 2 2 

^sp — l^spl ) Rsp — I'J^spl • (35) 

When Vsp(— oo) = K;p(oo), Tgp and Rgp coincide with |6spP and |aspP, respectively. 

We shall define the transmission coefficient T and the refiection coefficient R for 
the elementary excitations. Assuming that an excitation comes from the left, we write 
the asymptotic form of the wave function of the excitation as 

^^"^=i^^'(x)) 



(36) 



where ki and ^2 are real, and they satisfy 



e = qk + ^-{- + 2). (37) 

Equation ()37|1 is the Bogoliubov excitation spectrum in a uniform system where the 
condensate has supercurrent proportional to q [22] ■ Since Eq. (jHTj) is a fourth order 
equation for k, one can solve it and obtain ki and ^2 analytically. However, since we 
focus on low energy regions where a specific tunneling behavior appears, we only write 
the approximate forms of ki and /c2 for e <C 1 here, 

h ^ — — 77, (38) 

l + q 8(l + g)4' ^ ^ 

(39) 



1-q 8(1 -g)4' 
The amplitudes Uk and Vk are given by 



Uk 



k'^ 1 

1 + - + g - (jk Uqx+Sf,n{xm /.qN 

N 2(5 -gfc) ' ^^'^ 
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Vk 



1 + 



fc2 



e + qk 



(41) 



(42) 



\ 2{e-qk) ^ 
The Wronskian defined as 

yields a relation between a' and c\ which defines the transmission and reflection 
coefficients. One can easily prove from the Bogoliubov equations that W is independent 
of X when ip^ and '0' have the same energy. By evaluating W {ip''* , ip^) , one obtains 

J|2 





\' + 


I^fc2 






2 + 1 







+ c 



(43) 



which expresses the conservation law of the energy flux jTSj. It is obvious that the 
reflection coefficient R and transmission coefficient T are defined from Eq. ()43|) as 



R 
T 



-k2{\Uk2 


\' + 






ki{\uk^\ 


2+1 


Vki\ 


2) + g 



J|2 



(44) 
(45) 



When there is no supercurrent, R and T are simply given by \a \ and |c | , respectively. 



3.2. Calculation of transmission coefficient 

The condensate wave function of Eq. (PT|) takes the same form as that for a gray soliton 
except for the constant Xq. Consequently, the solutions of the Bogoliubov equations are 
also the same as those for a gray soliton which was obtained in Ref. Substituting 
Eq. (PT|) into the Bogohubov equations, one can analytically obtain four particular 
solutions. They are 



Unix 



^^gi[(fc„+g)x+sgn(x)|-] _ 
h 



2e 



7(a;j — 1 sgn(xj 



X 



q 



Vn{x) = A^e^[{'="-'/)^-«gn{x)|]l ^ 



^ J 7(x) + isgn(x) 



(46) 



X 



k 

g + -^(i-9^-7(a;)^-£) 



1.3 



(47) 



where kn satisfies Eq. fl37|) . Approximate forms of k^ and /C4 for e are given by 



k. 



3,4 



+ 2iv/l-g2 + 



qe . (1 + 2g2)£2 



+ 1 



1 - r 4(1 - g2 



(g + q^)e^ 



Equations (jHH|) and (jH^ show that {ui{x) , Vi{x))^ and (^2(3^), '^^2(3;))* describe scattering 
components. It is noted that there exist the scattering components whose energy is lower 
than the condensate potential, because the excitation spectrum is gapless in contrast 
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to the case of the superconductor. Equation (|iH|) shows that {us{x) , v^{x))^ on the left 
hand side of the barrier and {u4^{x), ^4(3;))* on the right hand side describe the locahzed 
components around the potential barrier, and {u4{x),V4{x)Y on the left hand side of 
the barrier and {u3{x) , V3{x))* on the right hand side describe the divergent component 
far from the potential barrier. The normalization constant A„ is expressed as 

^i^^ n = 1,3,4 



n 



where 



e " = 



4e + 2eqK + 2(1 - q^)kl + qkl 



AeJl + ^ + e - qkn 



+ i sgn(a;) 



^^r^kn{2e - 2qkn + kl) 



As 



1 + f +e 



qkn 



(49) 



(50) 



The normalization constant of the scattering components is determined to satisfy Eq. 
(ESI)- 

Two independent eigenf unctions of Eq. corresponding to two types of scattering 
process are obtained by omitting divergent components. One is the process in which an 
excitation comes from the left-hand side {ip\x)), and the other from the right-hand side 
{ip'^{x)). Here we consider the former written as 

The coefficients a', b\ c\ and d^ are the amplitudes of the reflected, the left localized, 
the transmitted, and the right localized components, respectively. They are functions 
of the energy e, the potential strength Vq and the phase difference y9. The boundary 
conditions at x = yield four equations to determine all the coefficients 

^'(+0) = ^'(-0), (52) 



^\x) 



X < 0, 
X > 0, 



(51) 



dilj^ 
dx 



+0 



d^^ 
dx 



+ 2W'(0.) 



(53) 



These equations are linear simultaneous equations for the coefficients a\ h\ c' and d\ 
and one can analytically solve them. Assuming e <^ 1 and Vq ^ 1, we can obtain 
analytical solutions of them within the first order of e or 

(3cosv9 + sin(y9 + A)e 



cosy? + 



— 4cos^(y9 + sin/? cos(/? — 2cos(y5 + 2 



2K) 



i(cosv9 -|- siny9)£ 



/ 3cos^V9 + 2cos(/9 — 1\ 
1 cosy? H — ) + 2smy9 e 



d' 



COS(y9 + 



4cOS^(y9 — Siny9 COS*/? + 2cOS(y9 — 2 



1 COSV9 E 



/(2^'), 



(3cos(j9 -|- sinoj -|- 4)5 ./ 3cos^o9 + 2coso9 — 1 
y.e + + 1 cos^ + 



(54) 
/{2z^), (55) 
(56) 
(57) 
(58) 
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The relations between the coefficients in ip'- and those in are 



When if is not equal to ipc-, we easily obtain the transmission coefficient T = |c' 
from Eqs. (jKHjl and It is expressed as 

where 



i|2 



^ = -W^2^ (60) 



cos^if 3cos'^v9 — 2cos<^ 

^ = V^ + — ^ — ■ 

The transmission coefficient has a peak at e = 0, and the peak has Lorentzian shape with 
half width A. This means that the potential barrier is transparent for the low-energy 
excitations. This behavior of the transmission coefficient, called anomalous tunneling, 
has been predicted for a current-free condensate IT^ . We can see from Eq. 
that the peak width decreases as (f approaches — f , and it becomes infinitesimal for 
{p ^ ipc- Thus, the anomalous tunneling is suppressed by the supercurrent. When ip 
reaches ipc, the transmission coefficient is expressed as 

^ ^ 4- (62) 

Obviously, the anomalous tunneling behavior does not exist any longer, and the 
transmission coefficient has only a small residual value. 

We plot the transmission coefficient as a function of e at ip = 0,|,y9c(— f),x 
and 71 are shown in Fig. 2, where we set Vq = 10. In Fig. IA2r a) we can clearly 
see the properties of the transmission coefficient mentioned above. In the region far 
from the anomalous tunneling peak, T increases with e in a conventional way and its ip- 
dependence disappears as shown in Fig. IA2r b). This is because the condensate potential 
terms in the Bogoliubov equations, written as I^E^oP? ^o^' begin to become so 

small compared to the kinetic energy term that the Bogoliubov excitations behave as 
single particle excitation. 

While we adopted the (5-function potential barrier of Eq. in the above 

calculations, the width of a potential barrier is finite in real systems. However, 
the treatment of the problem using the 5-function potential barrier gives us valuable 
qualitative insights, because it enables us to analytically calculate physical quantities 
such as the transmission coefficient. In fact, the result of the calculation of the 
transmission coefficient for the 5-function potential barrier ^3] is qualitatively the same 
as that for the rectangular potential barrier ^3] in the case of current-free condensate. 
Moreover, comparing Eq. ()28|) to Eq. ()A.16|) in Appendix, we can see that the relation 
between the supercurrent and the phase difference for the 5-function potential is also 
qualitatively the same as that for the rectangular potential. Accordingly, we consider 
that the (^-dependence of the transmission coefficient is valid not only for the (5-function 
potential barrier but also for a potential barrier with finite width. 
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4. Discussion 

Origin of the anomalous tunneling 

Since the anomalous tunneling exhibits behaviors specific to the resonant tunneling, it 
is expected that its origin can be attributed to the appearance of quasi-bound states 
with lifetime Xqb ~ Kagan et al. actually insisted that the quasi-bound states were 
induced by spatial changes of the condensate density. Their argument is as follows. The 
condensate density acts as a kind of potential for excitations. It inevitably reduces near 
the potential barrier, and consequently creates a potential well for excitations as shown 
in Fig. lAll Using an analogy from the resonant tunneling of single particles, such a 
potential well induces the quasi-bound state, and this is the origin of the anomalous 
tunneling. 

Here we try to elucidate mechanisms of the anomalous tunneling from another 
viewpoint of the prominence of the localized components with imaginary momenta of 
Eq. ()48|1 . Calculation of the probability density of each component in the left-incident 
state ip'' yields useful information to such viewpoints. Since the normalization condition 
for {u{x),v{x)y is given by 3\ 
1 
L 

the probability density of each component is defined as |ui(a;)p — for the incident 
component, |a'p(|M2(a;)p — |t'2(x)p) for the refiected component, |6'p(|'U3(x)p — |f3(x)p) 
for the left localized component, |c'P(|mi(x)P — |fi(a;)p) for the transmitted component 
or |c?'p(|-U4(x)p — |f4(x)p) for the right localized component. When Vq ^ 1 and e -C 1, 
one obtains the probability density at x = of each component, 

K(-0)p-K(-0)p~l + i^, (64) 
|a'P(M-0)p-M-0)P)~|af(i-i^), (65) 



dx{\u{x)\' - \v{x)\') = 1, (63) 



Ll\2 



3 + 3 cosip + sirnp\ 



bf{\us{-0)\'' - \vs{-0)[') - ^( - 2 - -^'^7 ^ (66) 



|c'p(K(+0)r-h(+0)r)~|c'r(i + i^), (67) 

mM+o)\^ - i..(+o)r) ~ ^( - 2 - ^_±i^^s^^±^y ^,s) 

We show the probability densities at x = of the transmitted component (TC) and 
the localized component (LC) in Fig. IA31 We can see from Eq. ()67|) and Fig. IA3r a) that 
the maximum value at £ = of the probability density of the transmitted component 
increases slightly as the phase difference increases. However, the maximum of the 
transmission coefficient does not change, because the probability density of the incident 
component also increases slightly as seen in Eq. When ip is not equal to ipc, the 

localized components prominently appear around e = where the anomalous tunneling 
occurs. When ip is equal to ipc, the localized components do not appear. Thus, the 
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localized components are induced only when the anomalous tunneling is effective. This 
result suggests that the localized components are crucial to understand the anomalous 
tunneling. It is to be noted that the probability densities of the localized components are 
negative Due to their negativity, the probability density of scattering components 
can be comparably large even near the potential barrier without totally raising the 
probability density of ip''{0). Accordingly, the localized components have a role to 
spread the scattering components across the potential barrier, and the appearance of 
the localized component is one of the origins of the anomalous tunneling. 



4-2. Observation of the anomalous tunneling 

In order to describe realistic parameters properly, we add the bars to the dimensionless 
parameters again as in the Eq. Q. Since the width of the potential barrier is finite in 
real systems, we consider a rectangular potential barrier, 

V{x) = Ve(^ - (69) 

where V and d are the height and the width of the potential barrier, respectively. The 
anomalous tunneling is predicted also for this rectangular potential barrier ^]. When 



the barrier is so high that e -C 1 is satisfied, the peak width of the transmission 
coefficient can be expressed as 



where 



(70) 



«^o = \ T^(V^-/i)- (71) 




One needs to create excitations with energy comparable to Arec for the observation of 
the anomalous tunneling in experiments. For that reason, the values of two parameters 
with length dimension are restricted. One is the barrier width d, which corresponds 
to thickness of the laser beam. The other is the system size, which is the axial length 
L of the system. The healing length ^ and the barrier intensity determine the 
restrictions. 

We consider the experimental setup of Ref. p/Ej. The healing length is ^ ~ l/xm 
in that experiment, because the total number of ^^Rb atoms, the system size and the 
radial harmonic trap frequency are ~ 3000, L ~ 80/im and 2t^uj± ~ 40kHz. We see 
from Eq. ()70|1 that the peak width decreases as the barrier width increases. This means 
that the potential barrier should be as narrow as possible. Since one can narrow down 
the laser beam waist, which corresponds to the barrier width, to the value comparable 
to the wavelength, the laser width can be ci ~ 1/im experimentally with barriers created 
by a blue detuned laser. 

We shall give the restriction on the system size. The anomalous tunneling occurs 
in a low energy region where the Bogoliubov spectrum is phononlike. Phonons with 
energy comparable to or smaller than Arec are available if the axial size of the system 
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L is larger than the wavelength Aa, where A a is easily calculated from Eq. (fTOIl as 
^ ~ KQ^e'^"'^. We set c? = 1 /im and ^ = 1 /im. With a barrier height of kq^ = 3, L 
should be comparable to or larger than 60 /xm. Since the system size is L ~ 70/im in 
the experimental setup of Ref. |T7], they barely overcomes this restriction. Moreover, 
one can enlarge the system size by changing the position of the optical end caps. Thus, 
the anomalous tunneling is expected to occur in the uniform BEC if one can use a 
sufficiently narrow laser beam to create the potential barrier. 

A procedure to observe the anomalous tunneling is as follows. We consider two 
BECs separated by a potential barrier. One can control the phase difference by moving 
the potential barrier at the velocity — where q is given by Eq. ()A.16jl (see Appendix). 
At first, one stimulates the condensate on the left hand side into phonon excitations with 
energy comparable or smaller than A^-ec by using the Bragg pulse. The phonon with 
wavelength A takes ^ to pass through the barrier, where Cg = is the Bogoliubov 

sound speed. After longer time than we detect the number of transmitted excitations 
Ut and that of refiected excitations Uj- by the time of fiight absorption images, and the 
transmission coefficient is given by T = nt/{nt + ^r)- Thus, the observation of the 
anomalous tunneling may be possible in an optimized setup. 

The dramatic (/^-dependence of the transmission coefficient discussed in Sec. Ill 
enables us to determine or estimate the phase difference between two condensates 
by measuring the transmission coefficient. The phase difference is near wk if most 
excitations at low energy transmit across the potential barrier, while it is near {n + |)7r 
if most excitations at low energy are reflected by the potential barrier, where n is an 
arbitrary integer. 

5. CONCLUSION 

In summary, we have investigated (^-dependence of the tunneling of the elementary 
excitations in a bosonic S-I-S junction by solving Bogoliubov equations, and we have 
found the signiflcant (^-dependence of the transmission coefficient. 

Calculating the probability density of the localized components, we have discussed 
mechanisms of the anomalous tunneling. We have shown that the localized component 
arises only when the anomalous tunneling appears; this means that the prominence of 
the localized components is one of the origins of the anomalous tunneling. 

We have discussed the feasibility of observing the anomalous tunneling in 
experiments. We expect that the anomalous tunneling may be observed in BECs trapped 
in a box-shaped potential. 
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Appendix: Josephson relation for a rectangular potential barrier 

We have already obtained the relation between the supercurrent of the condensate and 
the phase difference for a strong 5-function potential barrier given by Eq. ()28j) . and 
we have shown that the relation satisfies the Josephson relation. In this appendix, we 
shall calculate the relation between the supercurrent of the condensate and the phase 
difference for a rectangular potential barrier, and show that the Josephson relation is 
satisfied also in this case. 

We consider a BEC in the box-shaped trap mentioned in Sec. II. Instead of a 
(5-function potential barrier, we adopt a rectangular potential barrier with width d 
expressed in the dimensionless form as 

Vix) = Vei^-\x\). (A.l) 

We shall solve Eqs. ()12|) and (fT^ for this potential barrier and obtain the relation 
between q and (p. 

Outside the barrier, one obtains the solution of the same form as Eqs. (fTHj) and 

m, 

A' = ^{xf + q', (A.2) 

d r Q. 



S{x)-S{±-)= £jxj, 



T ^) + sgn(x) 



(A.3) 



where 



7(x) = ^l-qHanhUl-q'^{\x\-^ + Xo) \ . (A.4) 



Under the barrier, the general solution takes the form 

sn^{y/P^x, a) 
cn^{^/]3^x, a) 



^ - a H „; 7^1 (A.5) 



where 



a = A{0), (A.6) 



2 



= (A.7) 



P± = ° ^ ° ^, «:o = V2(V-/i). (A.8) 
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Here cn{u,q) and sn{u,q) are the Jacobi elliptic functions. Assuming KQ,e'^°'^ ^ 1, one 
can neglect the last term of Eq. (fT2|) . In this case, one finds the solution 



A' 



S 



\KoaJ 



■cosh2Koa; + 



- W 

\KoaJ 



(A.9) 



S 



d 

dx 



~ TT — 2 tan" 



(A. 10) 



The constants xq and a can be determined by the boundary conditions at x = ^, 



d 



'd 



d^r 



dx 



d^o 



j+o dx 



Substituting Eqs. (TOIl and (|Xl0|l into Eq. we obtain 

.2^ r /^(±|) 



<y9 = 7r — 2 tan 



-1 / K,oa 



+ grf + 2 



tan ^ 



7W 



tan 



(A.ll) 
(A.12) 

(A.13) 



Expanding Eqs. (jA.lljl . ()A.12|1 and ()A.13jl into power series of Kq^ and e one 
obtains 



2(1 + cosip) 

-e 2 , 



2g-Kod 
g ~ siiup. 



(A.14) 

(A.15) 
(A.16) 



We can clearly see from Eq. ()A.16|) that the relation between the supercurrent and the 
phase difference for the rectangular potential barrier takes the form of the Josephson 
relation as well as the case of the 5-function potential barrier. 
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Figure Al. Schematic picture of the bosonic S-I-S junction. 




Figure A2. The transmission coefficient T as a function of e and (/?, where Vq = 10. 
(a) Ts at low energy with ip — 0, ip — 7r/4, ip ~ ipc — 7i'/2, f = 37r/4 and (f = n are 
shown, (b) Ts are shown up to a higher energy region. 
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(b) 0.1 




Figure A3. The probability densities at a; = of the the transmitted component 
(TC) and the right locahzed component (LC) as a function of e and ip, where Vq = 10. 
(a) They are shown in a low energy region, (b) They are shown up to a higher energy 
region. 



